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The successful emulation of the Hubbard model 
[l| in optical lattices 0, has stimulated world 
wide efforts to extend their scope to also cap- 
ture more complex, incompletely understood sce- 
narios of many-body physics [H, 01 • Unfortu- 
nately, for bosons, Feynmans fundamental "no- 
node" theorem [6[ under very general circum- 
stances predicts a positive definite ground state 
wave function with limited relevance for many- 
body systems of interest. A promising way 
around Feynmans statement is to consider higher 
bands in optical lattices with more than one di- 
mension [7], where the orbital degree of freedom 
with its intrinsic anisotropy due to multiple or- 
bital orientations gives rise to a structural diver- 
sity, highly relevant, for example, in the area of 
strongly correlated electronic matter. In homoge- 
neous two-dimensional optical lattices, lifetimes 
of excited bands on the order of a hundred mil- 
liseconds are possible [H, but the tunneling dy- 
namics appe ars not to support cross-dimensional 
coherence [8l Il0|]. Here we report the first obser- 
vation of a superfluid in the P-band of a bipartite 
optical square lattice with S*-orbits and P-orbits 
arranged in a chequerboard pattern. This per- 
mits us to establish full cross-dimensional coher- 
ence with a life-time of several ten milliseconds. 
Depending on a small adjustable anisotropy of the 
lattice, we can realize real-valued striped super- 
fluid order parameters with different orientations 
P x ±P y or a complex- valued P x ± iP y order param- 
eter, which breaks time reversal symmetry and 
resembles the 7r-flux model proposed in the con- 
text of high temperature superconductors 
Our experiment opens up the realms of orbital 
superfluids to investigations with optical lattice 
models. 

Orbital optical lattices permit to prepare exciting new 
many-body scenarios with relevant counterparts in the 
area of strongly correlated electronic matter, where the 
orbital degree of freedom with its rich structure of possi- 
ble orientations gives rise to a wealth of interesting phe- 
nomena. For example, in the transition-metal oxides or- 
bital physics is believed to be a key element for under- 
standing their metal-insulator transitions, superconduc- 
tivity, or colossal magnetoresistance. Recent theoretical 
proposals have considered, for example, the formation of 
multifiavor and multiorbital systems d, jlj , Il3 1 . super- 
solid quantum phases in cubic lattices |14l . [la] , quantum 
stripe ordering in triangular lattices [l6j. or Wigner crys- 
tallization in honeycomb lattices (l7j . The recent predic- 



tion, that the life time of atoms in higher Bloch bands 
can be unexpectedly long, has further strengthened the 
potential significance of orbital physics in optical lattices 
[H, Previous experiments have demonstrated that 
the excitation of higher Bloch bands is in fact possible 
[Still- I n a recent experiment, investigating a homoge- 
neous quasi one-dimensional (ID) lattice of bosons, tran- 
sient partial coherence in the P-band has been observed 
[h| . In an extension to a two-dimensional square lat- 
tice, however, no cross-coherence could be established, 
which had previously been predicted in Ref. || as a 
consequence of the small transverse tunneling rate of the 
P-orbitals. This finding raises the question whether the 
implementation of non-trivial coherent orbital physics in 
the P-band, involving the freedom of orbital orientation, 
is at all possible in a homogeneous lattice potential. 

In this article we show that orbital superfluidity can 
be accomplished by an appropriate management of the 
tunneling junctions in a bipartite square optical lattice. 
We report the first observation of a 2D orbital super- 
fluid in the P-band. A key feature responsible for the 
formation of full cross-dimensional coherence is the com- 
position of the lattice of two types of sites with different 
well depths arranged as the black and white domains of 
a chequerboard. While the deeper wells provide local P- 
orbits, which mainly support longitudinal tunneling, the 
shallow wells contribute isotropic S'-orbits, which act as 
cross-dimensional tunneling links. By means of a popula- 
tion swapping technique we can selectively excite a large 
fraction of the atoms into the lowest P-band. While ini- 
tially the entire second Brillouin zone is smoothly popu- 
lated and the sample does not show any coherence, after 
a few milliseconds a long-lived P-band superfluid arises, 
which exhibits full coherence in both lattice dimensions. 
Only after times on the order of several ten ms coher- 
ence is lost as the atoms undergo collisions repopulat- 
ing the ground state S'-band. Depending on the adjust- 
ment of a small anisotropy in our lattice potential, which 
breaks the four-fold rotational symmetry of the P-band, 
we observe striped superfluid orders with different orien- 
tations or, if isotropy of the P-band is established, a gen- 
uinely complex-valued superfluid order, which exhibits 
staggered orbital currents. This latter case breaks time 
reversal symmetry, indicated by the presence of a sheet 
of staggered vortices with a unit flux quantum per pla- 
quette. It realizes a variant of the 7r-flux model [ll| pro- 
posed in connection with cupratc superconductors and 
represents a striking example beyond Feynmans no-node 
theorem. 

The core of our experimental set-up is a (quasi 2D) 
square optical lattice, composed of two classes of (tube- 






Li 


O 




£ 


□ 



step 2 




© © © © © 



FIG. 1: Population of excited bands, a, The lattice 
comprises two classes of lattice sites denoted by A and B. 
The grey area denotes the Wigner Seitz unit cell of the A- 
sublattice. b, Experimental sequence used to populate ex- 
cited bands, c, The populations within the first four Brillouin 
zones (BZ) are plotted versus the final value Of adjusted in 
step 2 of the sequence in b. The insets show illustrations of 
the first, second and forth Brillouin zones (upper row) and 
their observed populations for values of Of within the grey 
bar and at the dashed lines (lower row). 



shaped) lattice sites (denoted as A and B) arranged as 
shown in Fig. 1(a). The accordant optical potential is 
formed by crossing two optical standing waves oriented 
along the x- and y-axes with linear polarizations parallel 
to the z-axis, produced in the two branches of a Michel- 
son interferometer (details are deferred to the methods 
section). This allows us to adjust the time phase differ- 
ence between the standing waves 9 to any desired value 
(20j . For 9 < it/ 2 the ,4-sites are more shallow than the 
S-sites and vice versa. 

Two different kinds of anisotropics can occur in the op- 
tical potential. The first type, which affects the structure 
of the energy minima in the P-band and the correspond- 
ing superfiuid order, results from the fact that due to 
imperfect reflection optics the intensity ratios between 
the two traveling waves forming each of the two stand- 
ing waves are generally unequal. In order to control this 
anisotropy, appropriate polarization optics in one of the 
standing waves lets us rotate the linear polarization of 
the incident beam out of the z-axis by a small angle a 
but maintain the polarization of the reflected beam par- 
allel to the z-axis. Tuning of a allows us to adjust the 
structure of the energy minima in the P-band, which dis- 
plays a four-fold rotational symmetry if a = aj so ~ 7r/5. 
For a < cvi so or a > cvi so this symmetry is broken. A 
second kind of anisotropy results from an imbalance of 
the intensities in both interferometer axes, which locally 
is unavoidable due to the use of finite sized Gaussian 
beams. If a ~ a- lso , this anisotropy does not notably af- 
fect the structure of the P-band energy minima, where 
condensation takes place, and thus it is irrelevant for the 
structure of the superfiuid order observed in our experi- 
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FIG. 2: Evolution of band populations and momentum 
distribution, a, Mapping of band populations for increasing 
hold times in the lattice, b, Corresponding momentum spec- 
tra. The red dashed circles in the 12 ms spectrum highlight 
the higher order Bragg peaks. The lattice axes are oriented 
in the x and y-directions. The case a = is shown. 



ment. For details we refer to the methods section. 

We may populate excited bands via the following two- 
step population swapping procedure related with tech- 
niques used in Ref. [2l| (cf. Fig. 1(b)). Initially a Bose- 
Einstcin condensate (BEC) of Rubidium ( 87 Rb) atoms 
is prepared and the lattice potential is activated with 
9 = 9i ps 0.38 7r, such that the well depth of the S-sites is 
about 4.5 times larger than that of the .A-sites (cf. meth- 
ods section). A ground state lattice is thus formed with 
the shallow ,4-wells practically not contributing to the 
tunneling dynamics. In the second step, 9 is rapidly in- 
creased (within 0.2 ms, which is shorter than the nearest 
neighbor tunneling time) to a final value 9 = 9f > n/2 
such that now the .A-wells are significantly deeper than 
the £>-wclls (cf. Fig. 1(b)). Depending on the chosen value 
of 9f 1 population is transferred to excited bands accord- 
ing to the graphs shown in Fig. 1(c). Expectedly excita- 
tions to higher bands are not observed for 9f < tt/2. In 
this domain the ^-independent background population 
in the excited bands is a consequence of imperfections 
in our adiabatic band mapping procedure and a back- 
ground of uncondensed atoms arising outside the image 
plane because of a technical pecularity of our BEC pro- 
duction [22j| . For suitable values of 9f > n/2 the second 
(lower P-band) and the forth (lowest P/-band) become 
the most populated bands. Notably, the third band (up- 
per P-band) population always remains small. 

For the observations reported here, the population ra- 
tio between the second and the first band is maximized 
(9f = 0.53 it). After a variable hold time, we map the 
populations in the different Brillouin zones or record the 
momentum distribution of the atoms. We first studied 
the case a = 0, when the four- fold rotational symmetry of 
the P-band is broken. In Fig. 2(a) the distribution of the 
atomic population among the different Brillouin zones is 



3 



shown for hold times between 0.1 ms and 200 ms. After 
0.1 ms the second Brillouin zone is nearly uniformly pop- 
ulated. As time evolves, sharp peaks arise at the two fi- 
nite quasi-momenta K(__ + ) and K( + _), developing max- 
imal contrast at hold times around 12 ms and persisting 
for times up to 40 ms. Here, K/± .±) = h(±hk, ±Kk) with 
k = 2tt/X and A denotes the wavelength of the lattice 
beams. This indicates that a collision aided condensa- 
tion process occurs. The momenta K(_ )+ ) and K( +) _) 
are located on the boundary between the first and the 
second Brillouin zone. As the hold time is increased 
beyond a few ten milliseconds the sharp peaks disap- 
pear and the first Brillouin zone is repopulated, which 
is most clearly seen for long hold times above hundred 
milliseconds. In the corresponding momentum spectra 
(cf. Fig. 2(b)) past 2 ms hold time pronounced maxima 
at the momenta K(_ +j and K( + \ arise accompanied 
by higher order Bragg peaks (highlighted for the 12 ms 
picture by dashed red circles in Fig. 2(b)), which indicates 
the emergence of a finite momentum superfluid. The ab- 
sence of a zero momentum peak clearly discriminates this 
superfluid from the conventional superfluid observed in 
ground state lattices. When we carefully tune a in order 
to adjust isotropy of the P-band (the value of a has to 
match ai SO to better than 6a ps ±7r/100 (cf. methods 
section)) a different superfluid order emerges. As shown 
in Fig. 3(a) we then observe four condensation momenta 
K(_ +), K( +i _), K( +i _|_), and K(_ _) in the correspond- 
ing momentum spectra (and similarly in the Brillouin 
zone mappings). As is shown in Fig. 3(b) we may invert 
the initial anisotropy adjusting a > a- iso + 5a and recover 
condensation at two momenta, however with an orienta- 
tion rotated by 7r/2 with respect to the observations for 
a < ai SO — 6a. Comparable condensation time-scales and 
condensate fractions are found independent of the value 
of a. 

In Fig.3(c-g) the 12 ms momentum spectrum of 
Fig. 2(b) (replotted in Fig. 3(c)) is analyzed more quan- 
titatively. In Fig. 3(d) a side view of this spectrum along 
the direction denoted by S3 in Fig. 3(c) is shown. In 
Fig. 3(c) a sectional view along the direction denoted by 
Si in Fig. 3(c) is shown. The solid grey line, matching well 
with the observed data (black circles), shows a fit to a 
sum of the (tight binding) momentum Wannier-functions 
of the S- and P-bands, to account for the incoherent frac- 
tions in both bands, and of two equal Gauss-functions 
separated by ^/2hk, which account for the condensate 
fraction. The widths of the 5-band Wannicr-function 
and the double Gauss-function as well as the individual 
amplitudes of the three functions are used as fit param- 
eters. The contributions of the P- and 5-band Wannier- 
functions and the double Gauss-function are shown by 
the dashed red, dotted blue and solid green curves, re- 
spectively. The sectional view in Fig. 3(f) is along the 
direction denoted by S2 in Fig. 3(c). The grey line re- 
peats the incoherent fraction obtained from the data in 
Fig. 3(c). In Fig. 3(g) a section along the direction de- 
noted S4 in Fig. 3(c) is shown. Two Gauss-functions are 
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FIG. 3: Tuning of anisotropy, condensate fractions, a, 

Momentum spectrum for an isotropic P-band 01 — Qi so . b, 
momentum spectrum for the case of a > Qfi so . c, The 12 ms 
graph of Fig. 2(b) is replotted. d, Side view of the momentum 
spectrum in c along the direction S3, e, Sectional view along 
the direction Si in c. f, Sectional view along the direction S2 
in c. g, Sectional view along the direction S4 in c. 



used to fit the data yielding the coherent and incoher- 
ent fractions given by the solid green and dotted blue 
curves. Figs.3(e-g) clearly show that a substantial part 
of the atoms is condensed. 

The nature of the observed superfluid order parame- 
ters can be understood extending arguments discussed in 
Refs.Jl, Q. In the tight binding limit the wave function 
near a lattice site resembles a harmonic oscillator state. 
Since the atoms reside in the P-band, at the deeper wi- 
shes the possible states are superpositions of local P x and 
P y orbitals, which in absence of collisional interactions 
are energetically degenerate, if a = aj S0 . The repul- 
sive interaction favors a maximum angular momentum 
at these sites with an orbit configuration P x ± iP y cor- 
responding to angular momentum quanta ±S per atom 
[H, Q . Since the population swapping sequence ends with 
the lowest energy S-state of the shallow 2?-wells tuned 
into resonance with the P-state of the deeper A- wells, 
an approximate local S'-orbit is to be expected near the 
23-sites. In order to maximize intersite hopping, the local 
phases of adjacent orbits should match at their tunnel- 
ing junction. This suggests that the ground state within 
the P-band is characterized in configuration space by 
the complex order parameter illustrated in Fig. 4(a). A 
2A x 2A section of the lattice is shown with the deep A- 
sites and the shallow S-sitcs indicated by large and small 
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FIG. 4: Nature of superfiuid order, a, Orbital ordering 
within the P-band for b, Quasi-momentum rep- 

resentation corresponding to a (BZ = Brillouin zone). The 
four black disks denote the contributing quasi-momenta. c, 
Orbital alignment, if the lattice potential exhibits a small 
anisotropy (a < ai so ). d, Quasi-momentum representation 
corresponding to c. The dark (light) grey areas in a and 
c show the unit cell of the lattice potential (order parame- 
ters), e, Schematic representation of the phase ordering in 
c. The black lines show the orientation of plane wave fronts 
(subscripted by their Miller indices with respect to the prim- 
itive vectors ei,e2) compatible with the imposed phase or- 
der, f, Positions of the Bragg resonances (small red disks) 
in momentum space corresponding to the wave fronts in e. 
g, momentum spectrum calculated from the Bloch function 
corresponding to a and b. h, Lattice potential for 6 = 0.56 it, 
a = and a well depth parameter Vo = 6.2E rec - i, Bloch- 
function </>k (+ _% (a;, y) for the potential in h corresponding to 
the momentum components K( + _) and K(_ i+ i in d. 



grey disks. The plotted numbers denote the local phases 
of the wave function. The superposition of Bloch func- 
tions </>k (+ _j {x, y) + i </>k c+ +) {x, y) for the two inequiva- 
lent quasi-momenta K( + _) and K( ++ ) (derived from a 
band calculation for the potential in Eq. (1) of the meth- 
ods section including the lowest nine bands for a = a- lso 
and 9 = 0.53 it) reproduces this orbital order. The two 
pairs of equivalent quasi-momenta ~K(±,±\ and K( T ±) 
(plotted in (b) with respect to the first and second Bril- 
louin zones) represent the four degenerate energy min- 
ima of the lowest P-band, where the system is expected 
to condense (cf. Fig. 5 of methods section). The order in 
Fig. 4(a) has remarkable properties. With respect to the 
.A-sites it exhibits a pattern of staggered orbital angular 
momenta. With respect to the £>-sites it exhibits vortical 
plaquette currents and thus resembles the 7r-fiux phase 



discussed in the context of high T c -supcrconductors [11 1 
and the sta gge red vortex superfiuid, recently discussed 
in Refs. 



If a is adjusted below ai so , a difference in the energies 
of the P x ± Py-orbits is introduced, which may exceed 
the collisional energy gain of the maximum angular mo- 
mentum orbits. A band calculation for the case a = 
(and 6 = 0.53 tt) confirms that the energies of the band 
minima at K( + _) and K(_ + ) arc lowered with respect 
to those at K( + + \ and K(_ n by a few thousandth of 
the atomic recoil energy E ICC = ti 2 k 2 /2m (m denotes the 
atomic mass) corresponding to a few percent of the en- 
ergy width of the band (cf. methods section). The system 
then condenses into the striped state sketched in Fig. 4(c), 
which exhibits zero local angular momentum. Instead of 
local (P x ± iP y )-orbits (P x — P y )-orbits arise at the A- 
sites. The condensation momenta K( + _) and K(_ + j are 
shown in Fig. 4(d). An equivalent state composed of the 
quasi-momenta K( + + ) and K(_ \ (corresponding to lo- 
cal (P r +P y )-orbits) arises, if a > ai so . Remarkably, both 
order parameters in Fig. 4(a) and (c) break the transla- 
tional symmetry of the lattice potential, as is seen from 
their different unit cells plotted in the figure. 

The order parameters shown in Fig. 4(a) and (c) read- 
ily explain the higher order Bragg peaks observed in the 
corresponding momentum spectra. In Fig. 4(e) the B- 
sublattice is replotted for the case of the striped super- 
fluid order (Fig. 4(c)) with the red and blue discs denoting 
sites with opposite local phase. The black lines show the 
wave fronts of possible plane waves that can propagate 
through the lattice without suffering destructive inter- 
ference due to the imposed phase order. These planes, 
which are the lattice planes comprising only S-sites with 
the same local phase, are denoted by their Miller indices. 
In Fig. 4(f) the expected Bragg resonances in momentum 
space are plotted in agreement with the observations in 
Fig. 2(b). Similar arguments apply for the case of the 
staggered angular momentum order in Fig. 4(a). For this 
case we show in Fig. 4(g) the momentum spectrum calcu- 
lated by Fourier transforming the Bloch-function corre- 
sponding to the condensation momenta in Fig. 4(b). 

For the anisotropic lattice the observed momentum 
spectrum of Fig. 3(c) uniquely determines the order pa- 
rameter of Fig. 4(c). The isotropic case, however, requires 
a more careful discussion. The coexistence of regions 
with real-valued striped order but different orientations 
could in principle lead to a similar momentum spectrum 
as that of Fig. 3(a). Note, however, that an uncontrolled 
local anisotropy, which could lead to fragmentation into 
striped phases, is not available in our experiment and 
would not be compatible with the extreme sensitivity 
regarding the adjustment of a required to observe the 
spectrum in Fig. 3 (a). Considerations similar to those 
discussed, for example, in Ref. [26| suggest that for glob- 
ally degenerate condensation points, as realized in our 
experiment for a = ai so , a coexistence scenario with its 
excess tunneling energy at the phase boundaries is ener- 
getically unfavorable for small repulsive interactions, al- 
though a stringent theoretical proof is yet an open prob- 
lem. We thus may claim high plausibility for the conclu- 
sion that the momentum spectrum observed in Fig. 3(a) 
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in fact originates from the order parameter of Fig. 4(a). 
A unique proof, however, remains to be an exciting chal- 
lenge for future theoretical and experimental efforts. 

It is instructive to observe the important role of the 
local S-orbits in the shallow /3-wells. For simplicity, we 
consider the striped state. As is seen from Fig. 4(c) the 
5-orbits are essential for establishing coherence perpen- 
dicular to the direction of alignment of the (P x — P y )- 
orbits. Transverse tunneling between P-orbits is typi- 
cally an order of magnitude weaker than axial tunnel- 
ing [1, |i| . Tranverse tunneling between A- wells is thus 
expected to occur predominantly via the /3-wells. This 
picture is well confirmed, if the actual lattice potential 
(cf. Fig. 4(h)) is used to evaluate the Bloch- function (cf. 
Fig.4(i)) corresponding to the condensation momenta in 
Fig. 4(d). For better visibility of the P-orbits we have 
used 9 = 0.56 it here. An analog consideration applies 
for the staggered angular momentum order in Fig. 4(a). 
In fact, if resonant tunneling between A- and B- wells is 
suppressed by adiabatically switching off the /3-wclls af- 
ter the superfluid has been established, a rapid loss of 
transverse coherence is observed in the momentum spec- 
tra. Furthermore, a significant reduction of the lifetime 
of atoms in the P-band is observed. The latter finding is 
understood by noting that collisions within the /3-wells, 
where the atoms are in their local ground state, cannot 
reduce the P-band population. If, however, the atoms 
are constrained to reside exclusively in the .A-wells, en- 
hanced collisional decay to the lowest band is expected. 
Long P-band life times in a bipartite lattice have been 
previously predicted in Ref. [22j ■ 

In summary, we have reported the observation of or- 
bital superfluidity in the P-band of a bipartite optical 
square lattice, which exhibits full cross-dimensional co- 
herence with a life-time of several ten milliseconds. Our 
observations are explained by finite momentum super- 
fluid order parameters comprising S'-orbits and aligned 
P-orbits arranged in a pattern of staggered local phases. 
If the P-band is adjusted to be isotropic, the emerg- 
ing complex- valued superfluid order breaks time reversal 
symmetry, indicated by vortical plaquctte currents. Our 
experiment paves the way for exploring unconventional 
supcrfluids in optical lattices with potential relevance for 
material systems. We note that preliminary observations 
also show superfluidity with complex- valued order in the 
D- and P-bands. 



Methods Summary 

Optical Potential We couple a laser beam at a wave- 
length A = 1064 nm and I/e 2 radius u> = 100 /^m to a 
Michclson interferometer with its two branches crossed 
at a right angle in order to obtain the optical potential 

Vn — 2z ' 2 > 

— 4~ e ™° 1 71 cos ( a ) + y sin ( a )) elkx + e ~ z e ~ lkx ] 

+ e w z ( e lky +ee-' lky ) | 2 , (I) 




FIG. 5: Anisotropy and P-band energy minima, a, En- 
ergy surface E(q\, 52) of P-band plotted within the first Bril- 
louin zone for a = cti B0 . The red and blue discs indicate the 
local minima. Discs indicating equivalent minima share the 
same color, b, The energy difference SE of two inequivalent 
local energy minima is plotted versus Sa and Sr\. 



where 77 accounts for a small difference in the powers di- 
rected to both interferometer branches (a typical value 
for 77 is 770 ~ 0.95), e ~ 0.81 accounts for the power 
reduction in the retro-reflected beams due to imperfect 
optics, and the angle a permits us to adjust the degree 
of anisotropy introduced if e ^ 1. x,y,z are the unit 
vectors in the respective directions. An isotropic P-band 
arises if cos (a) w e, which amounts to a = a; S o ~ tt/5. 
The four local energy minima arising in the P-band at 
the edge of the first Brillouin zone (cf. Fig. 5(a)) are then 
degenerate. Tuning of a and rj affects the energy minima 
of the P-band in a notedly different way. Remarkably, if 
ol = aiso is adjusted, the structure of the P-band energy 
minima is not at all modified by tuning 77. This is shown 
in Fig. 5(b), where the energy difference SE of the two 
inequivalent local energy minima of the P-band (red and 
blue discs in (a)) is plotted versus Sa = (a — ai SO )/ai so 
and Sn = (rj — ?7o ) / *?0 ■ Within the dashed (red) lines 
we observe the momentum spectrum of Fig. 3 (a) while on 
the left and right the spectra of Fig. 3(b) or (c) are found. 
The insensitivity of the P-band minima against changes 
of 77 for Sa = is a highly relevant feature, since locally 
an imbalance of the standing wave intensities is unavoid- 
able due to the use of finite sized Gaussian beams. The 
time-phase difference 9 is proportional to the path length 
difference of the two interferometer branches, which is ac- 
tively stabilized with the help of a weak additional laser 
beam to better then 7r/200. A change of 9 by ir/4 can be 
obtained in less than 0.2 ms. 

Loading of Lattice A Bose Einstein condensate of typ- 
ically 2 x 10 5 87 Rb-atoms (in the P = 2,m F = 2 state) 
is produced in a magnetic trap with standard procedures 
at the crossing point of the interferometer branches. The 
magnetic trap is adiabatically deformed in order to pre- 
pare a nearly isotropic condensate fraction with 12 fxm ra- 
dius. Within 80 ms the lattice beam intensity is ramped 
up with a cubic spline to Vo/E lec = 6.2 with 9 adjusted 
to 0.38 7T. This leads to well depths of the A- wells and 
/3-wclls of 2.5-E l re c an d H-4 E Iec , respectively. Excitation 
of the P-band is obtained by ramping 9 within 0.2 ms to 
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a value of 0.53 71", thus switching the .A-wells and S-wclls 
to 7.5E lec and 5.0 E lec , respectively. 

Detection procedures Momentum spectra are recorded 
by switching off the lattice beams in less than 1/xs and 
subsequent free expansion of the atoms during 30 ms be- 
fore an absorption image is recorded. For Brillouin zone 



mapping the lattice beam intensity is exponentially low- 
ered with a time constant of 430 /is. Subsequently, bal- 
listic expansion during 30 ms and absorption imaging is 
applied. The images shown in Figs. 1, 2 and 3 are aver- 
ages over five individual measurements. 
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